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ABSTRACT 

We present maps of 7 young massive molecular clumps within 5 target regions 
in the C 18 [J— 1—0) line emission, using the Nobeyama 45m telescope. These 
clumps, which are not associated with clusters, lie at distances between 0.7 to 
2.1 kpc. We find C 18 clumps with radii of 0.5-1.7pc, masses of 470-4200 M , 
and velocity widths of 1.4-3.3 km s _1 . All of the clumps are massive and approxi- 
mately in virial equilibrium, suggesting they will potentially form clusters. Three 
of our target regions are associated with H II regions ("CWHRs" from Clump 
with H II Regions), and the other two are without H II regions (CWOHRs). The 
C 18 clumps can be classified into two morphological types: CWHRs shape a fila- 
mentary or shell-like structure, CWOHRs are spherical. The two CWOHRs have 
systematic velocity gradients. Using the publicly released WISE database, Class 
I and Class II protostellar candidates were identified within the C 18 clumps. 
The fraction of the Class I candidates among all YSO candidates (Class I+Class 
n) is > 50% in CWHRs, and < 50 % in CWOHRs. We conclude that effects from 
the H II regions can be seen in (1) spatial distributions of the clumps: filamen- 
tary or shell-like structure running along the H II regions, (2) velocity structures 
of the clumps: large velocity dispersion along shells, and (3) small age spreads 
of YSOs. The small spread in age of the YSOs show that the presence of H II 
regions tend to trigger coeval cluster formation. 
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Subject headings: ISM: kinematics and dynamics - - ISM: molecules — ISM: 
structures — ISM: H II regions — stars: formation 



INTRODUCTION 



Most stars, in particular high-mass stars (> 8M Q ), form as stellar clusters ( jLada fc Lada 

20031 ). Stellar clusters form in dense and massive molecular clumps (size ~ 1 pc, mass ~ 100- 



1000 M^, density ~ 10 3 5 cm 3 ) with l arge velocity widths (2-4 k m s x ) (jPlume et al.l 11997 



Ridge et al.l l2003[ lLada fc Lada! 120031 : iHiguchi et aD bood . boioh. Typical medi an ages of 
stellar members in these clusters range between 1 and 3 Myr (jLada Sz Ladall2003l ). although 
the reason of this spread is unclear. Since dense gas in the cl uster forming clumps is si gnifi- 
cantly dispersed by the feedback from newly formed stars (e.g.. lHiguchi et aU2009U2010l ). the 
initial conditions and details of cluster formation are not understood yet. Therefore, it is nec- 
essary to study molecular clumps in the early stages of cluster formation. Recently, Infrared 
Dark Clouds (IRDCs) have been identi f ied as dark extinction features against the bright 
emission using MSX (lEgan et al.l Il998l : ISridharan et al.l 120051) . and GLI MPSE 8/xm data 
(|Peretto fc Fulleii2009h . ! RDCs were catalogued by lSimon et al.l ( 120061 ) and lPeretto &: Fulleii 
(2009 ). Their physical conditions are similar to tho se of cluster-forming clumps ( iRidge et al. 
20031 : Lada fc Lada 2003 : Higuchi et al. 2009 , 2010 ). which are the dense regions in molecular 
clouds. Most IRDCs are found to be associated with YSO candidates in Spitzer and Herschel 
obser vations Jchurchwell et al.lbood l2007t beharveng et al.lbood . bold : IZavagno et aDl2007l . 
2010). Thus, IRDCs are suggested to be precursors of clust er-forming clumps which still 



prese r ve initial conditions for the formati on of stellar clusters ( lEgan et al.lll998l ; ICarey et al. 



19981 . 120001 : iRathborne et al.l I2006L l2010l ) because the dense gas component is not yet dis- 



persed by the stellar feedback. 



Deharveng et al.l I2005L l2010l ). In 



Some IRDCs are detected near H II regions (e.g. 
such cases, cluster formation may be triggered by H II regions. Recently, "bubbles" , which 
are surrounding expanding H II regions have been dete c ted from Spitzer and He rschel 
Jchurchwell et al.ll2006[ I2OO7I beharveng et al.llioosl . bold : Izavagno et al.lbo07l . boioh . The 
UV radiation from massive stars excite Polycyclic Aromatic Hydrocarbons (PAHs) at the 
bubble rims and also heat dus t grains. The former can be readily observed at 8/zm by IRAC 
dChurchwell et al.ll2006l . 120071) and the latt er at longer wavelengths using MIPS on Spitzer 



( iChurchwell et al. 



2007; Watson et al. 



2008). Cold dense material has been observed around 



the borders of the H II regions (e.g.. iDeharveng et al.ll2005l . l2010l : IZavagno et al.ll2007l . l2010l ). 



How the cold dense mater ial can be accumulated by the expanding shells from the H II region 



has been discussed (e.g., IDeharveng et al.ll2005l l2010l : IZavagno et al.l 120071 12010| ). Most of 
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the relevant observations have been conducted with mm dust continuum emission. Consid- 
ering dynamical effects such as shock waves and expanding H II region shells should strongly 
influence surrounding molecular gas and change the kinematic structures significantly. It is a 
subject of great interest to investigate velocity structures in molecular gas near the boundary 
of H II regions in a systematic way. 



Lopez-Sepulcre et al.l (120101 ) presented maps of C 18 0(J=2-1) toward 48 massive molec- 



ular clumps and found that a b out h alf of their sample sources have velocity gradients of 
~ 2 kms -1 pc -1 . iRagan et al.l (120121 ) found similar velocity gradients in NH 3 lines toward 
some IRDCs. These studies do not focus on the external environments around the clumps. 
Understanding the nature and origin of velocity gradients may be crucial in understanding 
the initial condition of cluster formation. If a clump is associated with an H II region, the 
velocity structure of the clump will be influenced by shock waves and an expanding shell. 
Understanding how clumps are influenced by H II regions require comparison of velocity 
structures among clumps associated with and without H II regions. 

In this paper, we present high-resolution observations (effective beam size of 22") with 
the Nobeyama 45m telescope toward seven nearby (< 2.1 kpc) massive clumps (e.g., infrared 
dark) within five target regions. Three of the objects are associated with well-known H 
II regions (named "CWHR", i.e., "Clump with H II Region"), the others are without H 
II region ("CWOHR", i.e., "Clump without H II Region"). We used the C 18 O(J=l-0) 
molecular line emiss ion which traces dense clumps with density around 10 3-4 cm -3 regions 
(jHiguchi et al.ll2009l ). C 18 is an excellent tracer of the H 2 column density because in most 
of cases they are optically thin, and the relative abundance of CO and its isotopes are in 
most cases unaffected b y the prevailing changing in the gas studied (see the discussion in 



Mauersberger fc Henkell (119931 )). We analyzed WISE archival data available for all target 



regions, in order to identify YSO candidates which may be associated with clumps, and 
estimate their relative ages through its Class I fraction in each region. Here, we investigate 
the differences of spatial distributions, physical parameters, velocity structures, and activity 
of star formation (e.g., YSO distributions, relative ages) between CWHRs and CWOHRs. 



2. OBSERVATIONS 



2.1. Source selection 



Our targe ts are selected from clumps that are expected to be in early stages of clus 



ter fo rmation (jEgan et al.lll998l ; iBeuther et al.ll2002l ; ISridharan et al.ll2005l ; iRathborne et al. 
20061 ). Our criteria are that the clumps are massive (~ 100 M ), which estimated from 
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dust continuum emission, and nearby within 2.1 kpc, in order to resolve the cloud structures 
corresponding to the Jeans length (assuming gas temperature of 15 K) or finer with the 45m 
telescope at 110 GHz. We finally selected seven young massive clumps (e.g., infrared dark) 
within 5 target regions, as listed in Table [TJ The association with H II region or lack of such 
association is also indicated in Table [TJ 



2.2. Observation with the Nobeyama 45m Telescope 



We carried out C 18 0( J=l-0; 109.782182 GHz) observations using the 45m telescope at 
the Nobeyama Radio Observatory (NRO) between December 2007 and May 2008. We used 
the 25-BEam Array Receiver System (BEARS) for the front end receiver, w hich has 5x5 



beams separated by 41". 1 w ith respect to each other on the plane of the sky (ISunada et al. 



2000l ; lYamaguchi et al.ll2000i ). At 110 GHz, the FWHM main beam size was 15" and the main 
beam efficiency was 0.45. For the back end, we used 25 sets of 1024 chann el auto-correlator s 
(AC), each with 32 MHz bandwidth and 37.8 kHz frequency resolution (ISorai et al.l 120001 ) . 
The frequency resolution corresponds to a velocity resolution of 0.10 km s -1 at 110 GHz. 
We applied on-the-fly mapping technique (jSawada et al.l 120081 ) and covered areas of 6'x6' 
to 10'xlO' which corresponds to 3.8x3.8 to 6.1x6.1pc at a distance of 2.1 kpc. During 
the observations, the system temperatures in the double sideband mode ranged from 400 
to 500 K. We used emission-free regions near the observed sources as the off positions. The 
standard chopper wheel method was used to calibrate the intensity scale into T£, the antenna 
temperature corrected for the atmospheric attenuation. The accuracy of telescope pointing 
was checked every 1.5 hr by observing SiO maser sources near (within ~ 45°) the target 
regions, and was confirmed to be better than 3". We correct the difference in intensity 
scales among the 25 beams of the BEARS by using reference data of W51, which was 
obtained in observations with a single sideband filtered SIS receiver S100, and the acousto- 
optical sp ectrometer (AOS) a s the backend. We used a convolution scheme with a spheroidal 
function (jSawada et al.ll2008l ) to calculate the intensity at each grid point of the final map- 
cube (R.A., Decl., and vlsr) data with a spatial grid size of 7". 5, which is half of the beamsize; 
the final effective resolution is 22" and the effective integration time at each grid point is 
about 15 minutes. As a result, the typical rms noise level in channel maps for all of the 
target regions was 0.14 ± 0.03 K in 0.1 km s _1 . 

We also carried out NH 3 (1,1; 23.694495 GHz) and NH 3 (2,2; 23.722633 GHz) observa- 
tions as the backup for bad weather. At 24 GHz, the FWHM main beam size and the main 
beam efficiency were 80" and 0.8, respectively. For the frontend, we used the 20 GHz HEMT 
receiver (H20). The standard chopper wheel method was used to convert the received inten- 
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sity into T£. The system noise temperature was 100 to 500 K. For the backend, we used the 
AOSs, which have a velocity resolution of 0.5 km s _1 at 24 GHz. We observed the NH 3 (1,1) 
and (2,2) lines simultaneously by using two AOSs each with a bandwidth of 32 MHz. The 
integration time of each point was typically 180 seconds, resulting in a rms noise level in the 
channel maps of 0.05 K. The telescope pointing was checked every 3 hours by observing the 
SiO maser sources. The pointing accuracy was better than 15". Data reduction and analysis 
were made with the NEWSTAR package developed at the Nobeyama Radio Observatory 
(NRO). 

3. RESULTS 

3.1. Spatial distribution and velocity structure 

Figures [T] to [5] show the integrated intensity maps of the C 18 O(J=l-0) emission (con- 
tours) superposed on infrared color images in which we used the WISE 3.4, 4.6, and 12 
jum composite color images, and WISE 22 /zm color images. We defined the C 18 emission 
enclosed the 3cr contours of the individual map as C 18 clump. For IRDC 20081+2720, 
AFGL 333, and G 79.3+0.3, the clumps are associated with well-known H II regions (see de- 
scription of individual sources), while for IRDC 19410+2336 and MSXDC G 053.11+00.05, 
there is no evidence of the existence of H II region associated with the clumps 0. Hereafter, 
we named the clumps associated with H II regions as "CWHR" (Clump with H II Region), 
and "CWOHR" (Clump without H II Region). 

From the morphological structure of the C 18 clumps, we can classify C 18 clumps 
into two types: a) filamentary or shell-like structures (IRDC 20081+2720, AFGL 333, and 
G 79.3+0.3), b) spherical (IRDC 19410+2336, MSXDC G 053.11+00.05). Although the 
number of samples is small, the CWHRs tend to have the shell-like structures, while the 
CWOHRs have spherical structures from our targets. In order to investigate the detailed 
velocity structures of the clumps, i.e., presence or absence of velocity gradients, areas of large 
velocity dispersion, we produced 1st and 2nd moment maps (Figures [71 to fTTT) . From the 1st 
moment images, we identify the objects which have a velocity gradient (~ 2 kms -1 pc _1 ) 
within the CWOHRs: IRDC 19410+2336, and MSXDC G 053.11+00.05. For the CWHRs, 
the velocity gradient within the clumps can not be well identified. Furthermore, from the 
2nd moment images, we can find that large velocity dispersion (er ~ 0.8 kms' 1 ) appears 
along the filamentary structure for AFGL 333 and IRDC 20081+2720, which may indicate 



They might be associated with ultra compact H n region within the peak emission of C 18 
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the interactions between molecular gas and the expanding shell of the H II region. For G 
79.3+0.3, there are two H II regions (G79.307+0.277: south-east, G79.29+0.46: north-west), 
we can see shell-like structures in the channel maps (0.25-1.2 km s _1 in Figure [TBI) of the 
whole of the G 79.3+0.3 clump. 



3.2. Physical properties of C 18 clumps 

In this Section, we describe the derivation of the physical properties of the C 18 clumps 
(i.e., the area within the 3 a contours). A summary of the physical properties is given in 
Table [TJ To define the clump radius, we estimated the projected area of the clump on the 
sky plane, A and the observed radius was derived as R b s =(A/Tc) 1 / 2 . Then we corrected for 
the spatial resolution as, 

-Rclump = [-Robs — (#FWHIVl/2) 2 ] 1//2 , (1) 

where #fwhm of 22" the full width half maximum (FWHM) effective beam size. The velocity 
width of the clump (one-dimensional FWHM velocity widths), AV^. mmp was derived as follows. 
We used the FWHM of the Gaussian fit of spectra averaged over the clumps to derive the 
virial masses of the clumps (see below). We derived the C 18 O(J=l-0) line widths, AVo DS 
by fitting the C 18 O(J=l-0) spectra averaged over the clumps with a Gaussian function. 
We estimated intrinsic velocity widths of the clumps, AVd ump , by correcting for the velocity 

1 /2 

resolution of the instruments, AV^ pec , as AVci ump =(AV r J 3S — AK, pcc ) . The velocity widths, 
AVd uni p range from 1.4 to 3.3 km s -1 , suggesting considerable non-thermal contributions; 
the thermal velocity widths (one-dimensional velocity widths in FWHM) of C 18 molecules 
are 0. 12-0.24 km s -1 with Tk = 10-30 K. The H2 column density of the clumps, A% 2 was 
calculated as, 

^h 2 = 4.7 x 10 13 ( -^-] " ( ^) exp ( 5 - 27 



1.7 x 10- 7 ) \K J \T CX /K 

V0.457 \l-exp(-T)J VKkms-7 1 J KJ 

where X c i8 is the fractional abundance of C 18 relative to H 2 , T ex is the excitation temper- 
ature of the transition, r is the optical depth, and f T^dv is the total integrated intensity of 
the C 18 line emission. For X c is , we used 1.7 x 10 -7 by Frerking et al. (1982). We applied 
t = 0.5, and r/(l — exp(— r)) of 1. Hofner et al. (2000) revealed that most of massive star 
forming regions are optically thin in C 18 emission, so that we consider that the assumption 
is reasonable. The total mass under the local thermodynamic equilibrium (LTE) condition, 
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M c \ ump is calculated as, 

M clump = 40 



i.7xio- 7 y V2ioopcy \ wo" J 

where D is the distance to the object and -R c ium P is the radius of the clump. The rotation 
temperatures, T rot (2,2 ; 1,1) were derived from the NH 3 data using the same method pre- 
sented by Ho & Townes (1983). We applied the relation between the kinetic temperature 
and the rotation temperature, Tk ~ T rot as proposed by Danby et al. (1988). Under the 
LTE condition, we consider that the excitation temperature is comparable to the kinetic 
temperature (T ex ~ T K ). We used the kinetic temperature derived from NH 3 dataE as the 
excitation temperature of the C 18 gas. We also calcu lated the virial m asses of the clumps 
using (M VIR /M ) = 209(it; <:lum p/pc)(Ay c i uni p/km s'M ^Ikeda et^ lw^. 



We obtained radii of 0.5-1.7 pc, masses of 470-4200 M , and velocity widths in FWHM 
of 1.4-3.3 km s -1 (see Tabled]). We also calculated virial masses, and found that all observed 
C 18 clumps are in virial equilibrium within the error in estimate including the uncertainty of 
abundance ratio and kinetic temperatures. In order to discuss the kinematics of the clumps, 
we calculated effective Jeans lengths of 0.12-0.34 pc. The effective Jeans length, Aj was cal- 
culated as Aj = 0.19 pc (T K /10 K) 1//2 {n^ n /10 4 cm" 3 ) where T K is the kinetic temperature 



of the clump and nn 2 is the density of the clump (IStahler fc Pallall2005l ). The density of the 
clump, nn 2 was calculated as nn 2 = (3/47r)(M c i ump //imH 2 i? C mmp 3 ), where M cmmp is LTE mass 
of the clump, /zm H2 is mean mass of molecule, and -R cmm p is the radius of the clump. We also 
calculated the effective crossing time, 2i? cmmp /0 v (<r v = AV r c i ump /v / 81n2) of (0.9-3) xl0 6 yr, 
and the internal kinetic energy, £ kin =(l/2)M cmmp cr v 2 of (0.17-6.6) x 10 46 erg (see Table |2]). 
Comparing the physical parameters between CWHRs and CWOHRs, there are no signifi- 
cant differences of the physical parameters (radii of 0.7-1. 7 pc, velocity widths in FWHM 
of 1.9-3.3 km s" 1 , masses of 740-3200 M in CWOHRs, and 0.5-1.4 pc, 1.4-3.0 km s~\ 470- 
4200 M Q in CWHRs) and kinetic energies ((0.6-6) xlO 46 erg in CWOHRs, (0.2-7) xlO 46 erg 
in CWHRs) 

W e calculated the viria l ratio using our results, and found that they are gravitationally 



bound. lLada fc Ladal (120031 ) suggest that the Star Formation Efficiencies [SFE=M c i uster / (M cluster -|- 



M c \ nmp )} of the cluster-forming clumps range from approximately 10 to 30 %. The SFEs of a 



2 Because we could not estimate the excitation temperature for low signal-to-noise ratio of the NH3 data 
(IRDC 19410+2336 and G 79.3+0.3), we adopted 15 K which is the typical temperature in IRDCs regions. 
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cluste r increases to a maximum value of ~ 30-40 % due to gas dispersal by the stellar feed- 
back (jHiguchi et al.ll2009l |2010| ; lLadaj 120101) . If we assume max i mum SFEs of 30 %, cluster 



masses are expected from 100 to 1000 M & . iKauffmann fc Pillail ( 120101 ) suggest an empirical 



mass-size threshold for cluster formation including massiv e stars. Using the mass-siz e rela- 
tion; M L T E (#ciump) > 87OM ( J R clump /pc) L33 presented in kauffmann fc Pillail told ) @, all 
of the C 18 clumps exceed the limit shown in Figure El suggesting that they have a potential 
to form a cluster including massive stars. 



3.3. YSO candidates in our objects 



In order to investigate the existence of YSO candidates in our target regions, we used the 
information on stellar members in the point source catalogs of the Wide-field Infrared Survey 
Explorer (WISE; Wright et al. 2 010; Jarrett et al. 2011) for 10'xlO' areas. We adopted 
the source classification scheme in lKoenig et al.l (|2012l). Thev applied the established source 
classification scheme of Spitzer IRAC data (IGutermuth et al.l 120081 . 120091 ) to the WISE four 
bands (3.4, 4.6, 12, and 22 /im). We selected all sources with photometric uncertainty 
< 0.2 mag in four bands. To reduce the risk that differential extinction will confuse our 
classification of YSOs, we have adopted color-color diagrams of [3.4]-[4.6] vs. [4.6]-[12] and 
[4.6]-[22] to distinguish between Class I and Class II candidates. Figure [T7] shows source 
color-color diagrams of [3.4]-[4.6] and [4.6]-[12] made using WISE 3.4, 4.6, and 12 jum data. 
Class I candidates are identified if their colors match: [3.4]-[4.6] > 1.0 and [4.6]-[12] > 2.0. 
Class II (T Tauri star candidates) are slightly reddened objects, and were selected by the 
condition of [3.4]-[4.6] > 0.25 and [4.6]-[12] > 1.0. Figure M shows the WISE 3.4, 4.6, and 



22jxm color-color diagrams. Class II sources are re-classifie d as [4.6]— [22] > 4.0 ( IKoenig et al. 



20121 ). The rest are considered to be diskless sources as in IKoenig et al.l (120121 ) . We selected 



"Class I candidates", which are selected from both [3.4]-[4.6] vs. [4.6]-[12] diagram in Figure 
QH and [3.4]-[4.6] vs. [4.6]-[22] diagram in Figure | 



We detected 52 Class I and 28 Class II candidates superposed on enclosed 3 o contours of 
the C 18 clumps. We expect that contamination from non-YSO candidates, such as galaxies 
(Stern et al. 2005) is small (less than 10 non-YSOs per square degree) from the analysis in 
Allen et al. (2007). Furthermore, comparing the spatial distributions of YSO candidates with 
C 18 emissions in channel maps, most of YSO candidates are associated with C 18 emission. 
Thus we presume YSO candidates superposed on C 18 clumps as a result of star formation 
within the clumps. Stellar candidates with intermediate to high masses were identified as 



3 i? c i ump is the effective radius of the clump, MLTE^ciump) is the mass of the clump 
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the Class I sources w ith the magnitudes lower than 5 in 12 fim (e.g., IZavagno et al.l 120071 ). 
Zavagno et al.l (120071 ) identified the brightest sources with Spitzer 8.0 fim < 6 mag, which 
corresponds to star that can produce Photo Dissociation Regions (PDRs). Most of extended 
nebulae are probably local PDRs created by the radiation of the massive stars. In addition, 
these are not massive enough to form H II regions, they can heat the surroundi n g dus t 
and create local PDRs. Although the wavelength is different from IZavagno et al.l ( 120071 ) . 
the positions of the bright sources with 12 /im < 5 mag corresponds to H II regions or 
PDRs in our objects. From these empirical results, 30 % of detected Class I candidates 
can be candidates of intermediate to massive s tars. The number de nsity of YSO candidates 

Lada fc Ladal koO$ ) (e.g 



does not satisfy the definition of "cluster" in ILada fc Ladal (120031 ) (e.g., 35 or more stars 
whose space densities exceed lM pc~ 3 ). Consequently, we call these YSO distributions as 
"asso ciation" in th i s stud y. Note that the AGN and PAH candidates have been removed 
as in iKoenig et al.l ( 120121 ). although there is still uncertainty because individual candidates 



have not been identified by the SED. Furthermore, in order to check the existence of more 
evolved YSO candidates, we use the point source catalogs of the 2MASS. We identified 
YSO candidates from the 2MASS J, H, and K s data which cover the mapping area s of ou r 
dense clump da t a, usi ng th e color-colo r diagr ams and methods as in ICohen et al.l ( 1l98ll ). 
Bessell fc Brettl ( 119881 ) , and ICarpenterl (120001 ) . Figure HH] shows J-H versus H-K s color- 
color diagrams of point sources for the all target regions with 10'xlO' regions. The spatial 
distributions of the YSO candidates are shown in Figures 171 to [TT1 

There is a difference in the spatial distributions of Class I and Class II candidates between 
CWHRs and CWOHRs. For CWOHRs (IRDC 19410+2336 and MSXDC G 053.11+00.05), 
the Class I candidates are distributed around 6-9 a contours of C 18 emission (within ~ 
0.5 pc), although it depends on observations. While, the Class II candidates are distributed 
~ 2pc apart from the positions of strong C 18 emission, which is close to edge of 3 o 
contours. For CWHRs (AFGL 333, IRDC 20081+2720), Class I candidates are distributed 
in the shell-like structure within the clump, i.e., they seem to align along the expanding 
shell from the associated H II region. For the whole clump of G 79.3+0.3, there seem to 
be scattering of the distributions of Class II candidates compared with the distributions of 
Class I candidates. 

The fraction of Class I candidates among YSO candidates, -Rciassi=^Vciassi/ (-^ciassi+^ciassi) 
was calculated in Table |3j -Rciassi was found to be higher than 50 % (50-60 %) in CWHRs, 
and smaller than 50% (10-40%) in CWOHRs (see Figure [20]). If we add YSO candidates 
identified by the 2MASS into calculation of the fraction, this tendency does not change. In 
particular, t he fraction of Clas s I candidates among YSQ candidates is p roportional to the 
cluster age ( lEvans et al.l l2009t iNakamura et al.ll2011t IKoenig et al.ll2012l ). Comparing the 
results of 2MASS with WISE, CWOHRs tend to have YSO candidates identified by 2MASS. 
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This result is consistent with small fraction of Class I candidates among YSO candidates 
in CWOHRs. The higher fraction of Class I sources indicates that CWOHRs are relatively 
younger than CWOHRs. Furthermore, this trend implies that the age spread of YSOs in 
CWHRs is relatively smaller than those of CWOHRs. 



DISCUSSION 



4.1. Comparison between CWHR and CWOHR 

There are three clear signatures that distinguish CWHRs from CWOHRs, (1) spatial 
distributions: shell-like or spherical, (2) velocity structures: absence or presence of the 
distinct velocity gradient, large velocity dispersion along shells, and (3) age spreads of YSOs: 
small or large (see Table H]). Our results imply that H II regions strongly impact on the spatial 
distributions, velocity structures, and age spreads of YSOs. However, such an impact is not 
evident in the other physical proprieties of clumps (size, velocity width, mass, and kinetic 
energies) . 

First, the shell-like structure in CWHRs indicates an interaction between the clumps 
and expansion shell from the H II region. Comparing the images of dust continuum (also for 
our targets) and molecular lines (e.g., CO(l-O), 13 CO(1-0)), the mole cular clumps associated 



with the H II region tend to have shell-like or filamentary structures (IDeharveng et al.l 12005 
20101 : iBrand et al.ll201ll ). These morphological structures may have been shaped by an H 
II reg ion. Filaments will be unstable due to the int er actions and will be easy t o fragment 



(e.g., iTeixeira et al.l 12006k 120071 : iPineda et al.ll201ll ). iFukuda fc Hanawal ( 120001 ) presented 
that there are two factors that can trigger star formation by an H II region; dynamical 
compression d ue to the H II region and cha nge in the self-gravity. In a numerical study of 
star formation, iHosokawa fc Inutsukal (120061 ) suggest that an expanding H II region is able to 
trigger star formation within the ambient molecular clouds traced by CO(1-0), 13 CO(1-0) 
lines. We discuss this in more detail in Sec J4.2.11 

Second, regarding the difference of velocity structures, shock waves and the expanding 
shell of an H II region should strongly influence s urrounding molecu lar gas and change 
significantly the kinematic structure of the clumps. IBrand et al.l ( 1201 ll ) present maps with 
various molecular lines around the H II region, Sh 2-217. They see no obvious velocity 
gradients. Velocity structures in CWHRs has not been discussed previously because most 
of previous surveys of dense material around H II regions have been conducted with mm- 
continuum. This is discussed in more detail in Sec J4.2.2l 



Third, the value of -Rciassi is considered to be a clear indication of the relative evolution- 
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ary s t atus of the YSO ass ociation (jEvans et al.ll2009l ; iNakamura et al.ll201ll ; iKoenig et al. 



20121 ). IKoenig et al.l ( 120121 ) present a plot of the ratio, .Rciassi f° r the cluster- forming regions 



as a function of the central cluster age. Using their correlation, the ages of our stellar asso- 
ciations obtained by WISE are less than 1 Myr. In particular, the ages of the associations 
are estimated to be less than 0.5 Myr if the value of -Rciassi is higher than 50 %. The higher 
fraction of Class I sources indicates that CWOHRs are relatively younger than CWOHRs. 
Furthermore, this trend implies that the age spread of YSOs in CWHRs is relatively smaller 
than those of CWOHRs. Comparing the spatial distributions of YSO candidates with the 
C 18 emissions, Class I candidates are distributed uniformly in the shell-like structure in 
CWHRs. On the other hand, the Class I candidates are distributed with the peak emission in 
C 18 Q, and Class II ca ndidates are distributed outward of the Class I candidates in CWOHRs. 



Zavagno et al.l ( 120071 ) discussed the triggering mechanism of an H II region, RCW 120 (d 



1.3 kpc), and identified YSO candidates using Spitzer data. The spatial distribution of YSO 
candidates within "condensation 1" around RCW 120, whose distance are ~ 0.8 pc from the 
exciting stars, are similar to those of G 79.3+0.3, IRDC 20081+2720, an d AFGL 333. Class I 
candidates are un iformly distributed and aligned on the shell structure (jZavagno et al.ll2007l ; 



Puga et al.l 120091). From th e spat ial distribution of the YSO candidates toward RCW 120 
(Figure 5 in I Zavagno et al.l (120071 )). we calculated the value of -Rciassi and found it to be ~ 
50 %, which is consistent with our results in CWHRs. 



4.2. Triggered cluster formation on the periphery of H II regions 

4-2.1. Triggering mechanism by Hn regions 



Star formation may be triggered by the expanding shell motion of an H II region for 
AFGL 333, IRDC 20081+2720, and the whole clump of G 79.3+0.3. Dynamical effects of H 
II regions can be seen in the spatial distributions, velocity structures, and small age spreads 
of YSOs. Here, we discuss a possible triggering mechanism, which leads these different 
signatures. 

There are mainly two models proposed for trig g ered star formation by H 1 1 reg ions, 



radiation-driven implosion (iLefloch &: Lazaret! 



collec t-and-collapse model ( Elmegreen fc Lada 19771 ; Zavagno et al 



1994 



Kessel-Devnet fc Burkertl 120031) and 



2006; 



Deharveng et al. 



2010l ). In the collect-and-collapse model, swept-up molecular gas is accumulated to form 
dense gas in the shock front by the supersonic expansion of an H II regi on. The shell is mainly 
occupied by dense molecular gas whose mass increases with evolution. JZavagno et al.l (12 007) 
discus sed the dynamical expansion of RCW 120 using the model of iHosokawa fc Inutsuka 
(120061 ). Figure 14 shows the gas shell swept up by the shock front, with densities in the 
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range of 10 cm . This layer may b ecome unstable and fragm e nt, which make dense co re 
(size ~ 0.1 pc, density ~ 10 4_5 cm~ 3 ; Elmegreen fc Ladal (ll9771 ); IWhitworth et al.l (119941 ) ) . 



This process tends to sweep up amounts of gas and allows to form massive dense clumps 
(IBrand et al.ll201ll ). which is considered to form clusters including massive stars. On the 
other han d, radiation-driven implosion models are p redicted to form only low- to intermedi- 
ate stars (IDeharveng et al.ll2010l ; IBrand et al.ll201ll ). 



Our clumps are massive enough to form a cluster including massive stars (IKauffmann fc Pillai 



2010l ). In addition, CWHRs have a few candidates of intermediate to massive stars are dis- 



tributed around the H II region in G 79.3+0.3 and AFGL 333. From those results, we 
suggest that the collect-and-collapse model may be more adequate than the radiation-driven 
implosion. 



4-2.2. Kinematic energy injection from H II regions 

Let us assume that there is one 20 M & star, which is located at 0.5 pc (the separations 
between the positions of H II region and the strong C 18 emission peak are 0.3 to 1 pc) apart 
from the strong C 18 emission peak of the clumps. We estimate the internal kinetic energy, 
-Skin of the individual molecular clumps to be ~ 10 45_46 erg (see Table [2]). We cons ider the 



kineti c energy transferred to the clump through the stellar wind. Using the values of lAbbott 



(119821 ). we estimated the kinetic energy supply rate by the stellar wind from a 20 M & star 
to be ~ 10 35 ergs -1 . Here, we assume that the stellar wind is isotropically driven to the 
molecular clump. The kinetic energy of stellar wind during the crossing time (£"sw hi Niwa 
et al. 2009) is estimated to be ~ 10 48 erg, which is larger than internal kinetic energies of 
the clumps. Therefore our estimation suggests that H II regions have potential to change 
the kinematics of the clumps located on the border of the H II region. We can present that 
the energy input from the H II regions through the shocks of stellar wind are sufficient to 
drive the internal kinetic energy of the clump. 



4.3. Modes of cluster formation: coeval or non-coeval 



Triggere d star formation by an H II region can explain small age spreads of YSOs. 
According to IZavagno et al.l ( 120071 ) , the dynamical expansion of an H II region (excited by 
a star of ~ 20 M ) reaches a size of lpc in 0.2 Myr. An H II region will accumulate dense 
gas and derive the formation of filaments on the border of the H II region. Dynamical 
interactions between an H II region and a filament will occur. This filament will become 
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gravitationally unstable due to the interaction. Once dense filament forms, fragmentation 
will occur along the filament to form dense co r es, whose typical s eparation are the order 
of the Jeans length (ITeixeira et al.l 120061 . 120071 ; iPineda et al.l 1201 ll ) . This mechanism will 
coevally form dense c ores along the fila ment. These cores will immediately start collapsing 
in the free-fall time (jlkeda et al.l 120071 ). which explains the small age spread of YSOs in 
CWHRs. In fact, fragment components can be seen in the channel maps around the H II 
region in Figure HH to [T6j YSO candidates with the separations of the order of Jeans length 
are distributed along the filament in observed CWHRs. 

For CWOH Rs, the clumps have sphe rical structure. They are more stable than filamen- 



tary one (e.g.. Ilnutsuka &: Miyamal Il992l ). Nev ertheless, a spherical clump will eve ntually 



fragment, which initiates cluster formation (e.g.. lHartmannll2002l ; lBonnell et al.ll2003l ). After 
fragmentation, the local condensations wi ll become dist i nct peaks, i .e., dense cores, where 
stellar clusters will be eventually born. iHiguchi et al.l (120091 120101 ) classified 14 cluster- 
forming clumps into three evolutionary stages and found that spherical clumps in the early 
stage tend to have a single peak. This result may imply that the fragmentation rate in a 
spherical clump can be lower than in a filament. During the cluster formation, the cores 
within the clump s will form stars in a wide mass range and dense gas is dispersed by their 
stellar activities ( IHiguchi et al.l 120091 . 120101 ). Once clusters are formed, young stars, par- 
ticularly most ma ssive ones will disp e rse th e surrounding gas with evolution making cavi- 
ties around them (jFuente et al.lll998l . 120021 ) . Stellar feedback also forms dense cores (e.g., 



Knee fc Sand~el3 boooh . the structures of the evolved clumps are expected to be less con- 
densed than those of the younger clumps. The differences of the fragmentation rate between 
CWHRs and CWOHRs will explain the difference of the age spreads of YSOs. 



Pomares et al.l (120091 ) observed swept up dense gas surrounding Hn region, RCW 82 



and several Class I candidates along the expanding shell. iPuga et al.l (120091 ) investigated the 
spread of ages in Sh-284 with Spitzer/IRAC, and presented that despite the different scales 
(~ 20 pc) observed in its multi ple-bubble morph ology, the spread of ages among the powering 
high-mass clusters is ~ 3 Myr. iBik et al.l (120121 ) presented the evidence of an age spread of at 
least 2-3 Myr in W3 Main. Those scales of H II regions are a magnitude l arger than those o f 
our objects. If we use the age spread-linear size relationship described by lElmegreenl ( 120001 ). 
the age spread is roughly estimated to ~ 0.5 Myr (i.e., consistent with free-fall time) in our 
observed scales. This value is consistent with the ages of associations derived from relation 
between the ratio of Class I candidate and cluster age in iKoenig et al.l (120121 ) . Lifetimes of 
the low to intermediate stars (Class I and Class II candidates) are considered to be ~ 10 Myr 
( iStahler fc Pallal 120051 ). Thus, an age spread of 0.5 Myr is significantly small to speak of 
"coeval" cluster formation. 
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We propose that there are two modes of cluster formation; one is coeval cluster forma- 
tion, the other is non-coeval one. This difference is related to their environments (existence 
or non-existence of the H II region) of the associated clumps. H II regions play three roles: 
(1) they accumulate the dense gas by an expanding shell, (2) they shape the dense compo- 
nent into a filament, and (3) they inject energy. Our study can reveal differences between 
the mode of cluster formation as we focused on clumps in the earliest stages of cluster forma- 
tion. Interestingly, both CWHRs and CWOHRs have potential to form a cluster containing 
massive stars. In particular, CWOHRs contain candidates of intermediate to massive stars, 
although their fraction of Class II candidates is relatively higher. If the differences of cluster- 
forming modes depend on the differences of the environments of the clumps, some unknown 
issues of cluster formation (e.g., IMF, star formation efficiency, age spread of the cluster 
members) can be clarified by carrying out such an investigation toward the clumps in the 
early stages of the cluster formation (e.g., IRDCs). Although more statistical studies should 
be needed, these studies will lead the understanding the fundamental topic of the cluster 
formation. 



5. CONCLUSIONS 

We have mapped 7 young massive molecular clumps within 5 target regions in the C 18 
line emission with the Nobeyama 45m radio telescope. Our aim of this study is to understand 
the dynamical impact of H II regions on the clumps. We found different signatures of the 
clumps with H II region as compared to clumps or without one. Our results and conclusions 
are summarized as follows: 



1. We made the C 18 maps with a size of ~ 6'x6' to 10'xlO' areas of 7 nearby (D < 
2.1 kpc) young massive clumps within 5 targets. As a result, all the clumps are detected 
in C 18 lines, whose radii, masses, and velocity widths are 0.5-1.7pc, 470-4200 M Q , 
and 1.4-3.3 km s , respectively. They are similar to those of the cluster-forming 
clumps. 

2. All the clumps are likely to be in virial equi librium, and exceed the m assive star 



formation limit; mass-size relation presented in lKauffmann fc Pillail ( 120101 ) . suggesting 



that they have a potential to form a cluster including massive stars. 

3. In order to investigate the existence of YSO candidates in our target regions, we iden- 
tified stellar members using the point source catalogs of the WISE. There are 52 Class 
I and 28 Class II candidates within the all C 18 clumps. 30 % of the detected Class I 
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candidates may be intermediate to massive stars. The fraction of Class I to YSO can- 
didates: -Rciassi=N C iassi/ (A^ciassi+^ciassn) is presented in Table [3j It is > 50% (50-60 
%) in CWHRs, and < 50% (10-40 %) in CWOHRs. 

4. There are three clear signatures that distinguish CWHR and CWOHR, (1) spatial 
distributions: shell-like or spherical, (2) velocity structures: absence or presence of the 
distinct velocity gradient, large velocity dispersion along shells, and (3) the age spread 
of YSOs: small or large. 

5. An H II region can play three roles: (1) accumulate the dense gas by expanding shell, 
(2) shape the dense component into a filament, and (3) inject energy. The expanding 
shell can accumulate dense gas around the border of H II region and make the dense 
gas filament. Interaction between H II regions and filaments will tend to form dense 
core coevally. Once dense cores form, these cores will immediately start collapsing in 
a free-fall time, which can explain the small age spread of YSOs. 



We thank the anonymous referee for constructive comments that helped to improve this 
manuscript. We are grateful to the staff of the Nobeyama Radio Observatory (NRO)0 for 
both operating the Nobeyama 45m telescope and helping us with the data reduction. We 
also thank S.Komugi for useful discussion. 



A. Individual sources 

A.l. MSXDC G 053.11+00.05 

From the WISE composite image (Figured]), we can see two main YSO associations in 
MSXDC G 053.11+00.05 region located on north and south (also see Figure[7|). The channel 
maps in Figure [T21 clearly show that strong C 18 line emission is associated with these YSO 
candidates in different velocity ranges, for the south one which is brighter is in a range of 
20.8 km s _1 to 22.8 km s -1 , for the other in north is in a range of 20.8 km s _1 to 24.2 km 
s _1 . The C 18 clump seems to have a velocity gradient along the north-south direction. 
We consider that these emission components correspond to parental clumps of these YSO 
associations (see Figure [7]), whi ch are labeled as CI for the south one and C2 (corresponds 



to MSXDC G 053.25+00.04 in iRathborne et all (120061 )) for the north one. For clump CI 



4 Nobeyama Radio Observatory is a branch of the National Astronomical Observatory of Japan, National 
Institutes of Natural Sciences. 
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compact dust condensation is detected by iRathborne et al.l (120061 ) in 1.2 mm continuum 
observations. MSXDC G 053.11+00.05 CI have region of large velocity dispersion where 
YSO candidates are distributed, while MSXDC G 053.11+00.05 C2 have the large velocity 
dispersion where YSO candidates are not distributed. 



A. 2. IRDC 19410+2336 



From the images of Figure [21 we can see a YSO association in the WISE compos- 
i te image around IRDC 19410 +2336 central region, corresponding to IRAS 19410+2336 
(jMartm-Hernandez et al.l 120081 ) . The channel maps in Figure H3] clearly show that strong 
C 18 line emission is associated with YSO candidates in a velocity range from 21.2 kins" 1 
to 23.2 kms" 1 , and the C 18 clump seems to have a velocity gradient in the range from 20.8 
kms" 1 to 23.8 kms -1 in the direction from north-east to south-west. The multiple CO(J=2- 
1) outflows (blue; 5-18 kms " 1 , red; 26-47 km s" 1 ) have been detected in the direction from 
north-east to south-west by iBeuther et al.l (120031 ) . which implies that IRDC 19410+2336 is 
a young active star forming region. 



A.3. IRDC 20081+2720 



There is bright infrare d emission in the WI SE composite image (Figure [3]), which cor- 
responds to an H II region ( IBeuther et al.l 120021 1 . The H II region is located ~ 0.3 pc away 
from the strong C 18 line emission. In a velocity range from 5.2 kms -1 to 5.8 kms" 1 , we 
can see the C 18 emission have shell-like structure, which may be shaped by the H II region, 
i.e., from the interaction between dense gas and expanding shell. The YSO candidates are 
distributed within the shell-like structures of C 18 emission around the H II region. 



A.4. G 79.3+0.3 



There is a bright H II region (Figure HD, DR 15 (G 79.307+0.277), which lies behind 
and slightly south (~ 0.6 pc) of the G 79.3+0.3 region (IRedman et al.l 120031 ) . DR 15 is 
formed by one or two B0 zero-age main-sequence stars, observed as the far-infrared source 
(Odenwald et al. 1990), having a total luminosity of ~ 3 x 10 6 L & (Oka et al. 2001). A ring 
shaped nebula in 22 /j,m emission corresponds t o a photo-dissocia tion region, G79.29+0.46, 
which is the outskirts of an evolved massive star (IRizzo et al.ll2008l ). From the morphological 



distributions between C O and 22 /im emission, in a velocity range from —0.25 kms to 
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1.2 kms -1 , we can see the C 18 emission have an shell-like structure in the north area, 
which is an possibly compressed by G79.29+0.46. The YSO candidates are distributed in 
a shell-like fashion around the H II region (DR 15). The channel maps in Figure [TBI clearly 
show that strong C 18 line emission is associated with YSO candidates for the east one in a 
range of 0.25 kms -1 to 2.8 kms" 1 . For the other condensations in west, it is still dark in the 
22 fim image. In G 79.3+0.3 region, the global trend of velocity gradients can not be well 
defined. The c l ump o f the association for the east one are labeled as CI (corre spond to PI in 



Redman et all (120031 ) ). For the west one labeled as C2 (correspond to P3 in iRedman et al. 



( 20031 )). the C 18 clump is associated with an infrared dark area 



A. 5. AFGL 333 



We can see several YSO candidates in the C 18 emission (Figure [5j, whose separations 
are comparable to a Jeans length (~ 0.2 pc). In addition, we also see the bright infrared 
emission m the east of the C 18 e mission , which corresponds to an H II region generated by 
a B0. 5 type star (jHughes fc Vinerlll982l ). It is located ~ lpc away from the strong C 18 
line emission. In the north area of the C 18 emission, we can see t he relatively brigh t and 
diffuse 22 fivcv emission, which corresponds to a Herbig Ae/Be star (ISakai et all 120071 ) . The 
channel maps in Figure [TH] clearly show that strong C 18 line emission is associated with 
these YSOs candidates in a range of —50.2 kms -1 to —47.2 kms -1 . In a velocity range from 
—48.2 kms" 1 to —49.2 kms -1 , the C 18 emission have an shell-like structure. The YSO 
candidates are distributed in a shell-like fashion around the H II regions. In a velocity range 
from —50.8 kms -1 to —48.2 kms -1 , we can see the C 18 emission in northern edge have 
also an shell-like structure, which seem be compressed by the radiation or stellar wind of the 
Herbig Ae/Be star. 
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19 h 30 m 00" 29^45' 39 m 30" S9 m 15 B S9 m 00 s 38 m 45 s 19 h 30 m 00 s 29 m 45 s 29 m 30 s 29 m 15 5 29 m 00 s 2B m 4S s 



Right Ascension (J3000) 

Fig. 1. — [Left] Integrated intensity map of the C 18 0( J =1-0) emission (contours) superposed 
on the WISE 3.4 /mi, 4.6 /mi, and 12 /mi composite color image for MSXDC G 053.11+00.05. 
[Right] Integrated intensity map of the C 18 O(J=l-0) emission (contours) superposed on the 
WISE 22 /mi color image. The contours with the intervals of the 3 a levels start from the 
3 a levels for integrated intensity maps (la is 0.53 K km s _1 ). The pink line shows the lpc 
scale. The filled circle at the bottom left corner in each panel shows the effective resolution 
in FWHM of 22". 
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Right Ascension (J2000) 

Fig. 2.— Same as Figured] for IRDC 19410+2336. la is 0.26 K km s" 1 . 
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Fig. 3.— Same as Figured] for IRDC 20081+2720. la is 0.35 K km s _1 . 
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Fig. 5. — Same as Figured] for AFGL 333. la is 0.62 K km s . 
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Fig. 6. — Plot of r adius vs. mass and the mass-radius relation (solid line) derived in 
Kauffmann fc Pillail (120 lOh . The mass-radius relation up to which clumps form massive 
stars. 
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Right Ascension (J2000) 

Fig. 7. — (Left) The integrated intensity map of the C 18 O(J=l-0) emission (contours) 
superposed on the 1st moment map (color). (Right) The integrated intensity map of the 
C 18 O(J=l-0) superposed on the 2nd moment map (color) for MSXDC G 053.11+00.05. 
Small circles show Class I candidates, cross marks show Class II candidates, and large squares 
show YSO candidates identified by 2MASS. 
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Fig. 8.— Same as Figure for IRDC 19410+2336. 
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Fig. 9.— Same as Figure for IRDC 20081+2720. 
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Fig. 10.— Same as Figure for G 79.3+0.3. 
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Fig. 11 — Same as Figure Hfor AFGL 333. 
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Fig. 12. — The velocity channel maps of the C 18 O(J=l-0) emission with contours and 
grayscale for MSXDC G 053.11+00.05. The contours with the intervals of the 3 a levels 
start from the 3 a levels, where the 1 a noise levels are 0.14 K in T£. The velocity intervals 
are 0.4 km s _1 . 
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Fig. 13. — Same as Figure [T21 for IRDC 19410+2336. The contours with the intervals of the 
3cr levels start from the 3 a levels, where the 1 a noise levels are 0.11 K in T£. The velocity 
intervals are 0.4 km s" 1 . 
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Fig. 14. — Same as Figure fT2l for IRDC 20081+2720. The contours with the intervals of the 
3 a levels start from the 3 a levels, where the 1 a noise levels are 0.14 K in T£. The velocity 
intervals are 0.4 km s" 1 . 
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Fig. 15. — Same as Figure [T2l for G 79.3+0.3. The contours with the intervals of the 3 a 
levels start from the 3 a levels, where the 1 a noise levels are 0.20 K in T£. The velocity 
intervals are 0.4 km s _1 . 



02 28 45 30 15 00 30 

RIGHT ASCENSION (J 2000) 

Fig. 16. — Same as Figure [T2l for AFGL 333. The contours with the intervals of the 3 a levels 
start from the 3 a levels, where the 1 a noise levels are 0.12 K in T£. The velocity intervals 
are 0.4 km s _1 . 
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Fig. 17. — WISE 3.4 /xm, 4.6 /im, and 12 /im color-color diagram proposed by Koenig et al. 
(2012) for 10'xlO' areas. Right side of dashed lines indicate the boundaries between Class 
I and Class II candidates, and the left side of dashed lines indicate diskless candidates and 
Class II candidates. 
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Fig. 18. — WISE 3.4 /zm, 4.6 //m, and 22 /im color-color diagram proposed by Koenig et al. 
(2012) for 10'xlO' areas. Top of dashed lines indicate the boundaries between Class I and 
Class II candidates, and bottom lines indicate the boundaries between Class I and transition 
disk candidates. 
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Fig. 19. — J-i? vs. H-K s color-color diagrams of 2MASS point sources for our target objects 
for about 10'x 10' areas. The solid curves show the locus of un-reddened main sequence (pink) 
and giant stars (blue) in the CIT system (Bessell & Brett 1988). The photometric system 
of the intrinsic colors for Dwarfs (pink) and Giants (blue) corresponded to that of 2MASS 
system. 
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Fig. 20. — Bar graph of the fraction of Class I (pink) and Class II (yellow) candidates listed 
in Table |3j The number shows iV classI and iVciassii, which are enclosed within the 3 a emission 
of C 18 lines contours. 



Table 1. Physical parameters of the C 18 clumps 



Source Name 


RA (J2000) 


Dec (J2000) 


D [kpc] 


-RclumpM 


AV clump [km s- 1 ] 


m V ir[m s ] 


N(H 2 )[x 10 22 cm- 2 ] 


M clump [M ] 


Tex IK] 


H II region 


MSXDC G53. 11+00.05 CI 


19h29ml8s 


17d56m33s 


1.8 


1.2+0.4 


1.9+0.3 


920+560 


3.4+0.7 


3200+1400 


18 


No 


MSXDC G53. 11+00.05 C2 


19h29m37s 


18d02m00s 


1.9 


0.7+0.2 


2.0+0.3 


610+360 


2.4+0.6 


740+370 


15 


No 


IRDC 19410+2336 


19h43ml0s 


23d45m04s 


2.1 


1.7+0.5 


3.3+0.3 


4000+2000 


1.7+0.5 


2900+1200 


15 


No 


IRDC 20081+2720 


20hl0ml3s 


27d28ml8s 


0.7 


0.6+0.1 


1.4+0.3 


230+140 


2.6+0.4 


470+150 


15 a 


Yes 


G79.3+0.3 CI 


20h32m23s 


40d20m00s 


0.8 


0.7+0.1 


2.2+0.3 


650+310 


4.2+0.5 


1000+340 


15 a 


Yes 


G79.3+0.3 C2 


20h32m06s 


40dl8m00s 


0.8 


0.5+0.1 


2.8+0.3 


860+360 


3.9+0.5 


630+230 


15 a 


Yes 


AFGL 333 


02h28ml9s 


61d29m40s 


2.0 


1.4+0.4 


3.0+0.3 


2500+1300 


3.7+0.4 


4200+1800 


15 


Yes 



Note. — This table shows the name, center coordinates, distance, mass of the clusters, R c i ump - radius, AV c | ump : velocity width, iV(H2): mean H2 column density, -M clump : clump mass, T cx : 
excitation temperature derived by kinetic temperature of NH3 of the clumps. 

a Because wc could not estimate the excitation temperature for low signal-to-noise ratio of the NH3 data, we adopted 15 K which is the typical temperature in IRDCs regions. 



I 

CO 

00 

I 



-39 - 



Table 2. Kinetic energies of C 18 clumps 



Source Name 


[pc] 


Tcross [yr] 


Ekin [erg] 


Ekin/T"cro SS [erg S X ] 


MSXDC G53. 11+00.05 CI 


0.21 


3xl0 6 


2.1X10 46 


2.2xl0 32 


MSXDC G53.11+00.05 C2 


0.17 


2xl0 6 


5.6 xlO 45 


l.lxlO 32 


IRDC 19410+2336 


0.34 


2xl0 6 


5.9xl0 46 


7.9 xlO 32 


IRDC 20081+2720 


0.16 


2xl0 6 


1.7xl0 45 


2.9xl0 31 


G79. 3+0.3 CI 


0.13 


lxlO 6 


8.6 xlO 45 


2.0xl0 32 


G79. 3+0.3 C2 


0.12 


9xl0 5 


8.7xl0 45 


3.1 xlO 32 


AFGL 333 


0.20 


2xl0 6 


6.6 xlO 46 


9.8xl0 32 



Note. — This table shows the source name, the effective Jeans length of the clumps, 
the effective crossing time, the internal kinetic energy of the clumps, and the turbulence 
dissipation rate of the clumps. 



Table 3. Information of YSO candidates within enclosed 3 a contours (except within H II 

region) 



Source Name N classI N ciassn N 2M ASS N classI /(N classn +N classI ) [%] 



MSXDC G53.11+00.05 CI 4 (6) 3 (7) 3 57 

MSXDC G53.11+00.05 C2 3 (4) 4 (7) 43 

MSXDC G53.11+00.05 (total) 7 (10) 7 (14) 3 50 

IRDC 19410+2336 6 9 (11) 3 40 

IRDC 20081+2720 7 3 70 

G79.3+0.3 CI 11 9 55 

G79.3+0.3 C2 3 100 

G79.3+0.3 (total) 14 9 61 

AFGL 333 18 (21) 2 90 



Note. — Nciasaj: the number of Class I candidates within the clump, Nci ass n : the number of Class 
II candidates within the clump. The numbers in parentheses show total numbers of YSOs candidates, 
which arc distributed enclosed 3 a contours and distributed adjacent 3 a contours. N2MASS : the number 
of 2MASS YSO candidates within the clump, and Nc lassI /(NQ lassI +Nc lassII ): the fraction of Class I 
to the YSO candidates. 
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Table 4. Score sheet of different features between CWHRs and CWOHRs 



CWHRs (H II region) CWOHRs (Non-H II region) 



Source Name IRDC 20081+2720, G79.3+0.3, AFGL 333 MSXDC G53. 11+00.05, IRDC 19410+2336 

Spatial distribution filamentary, shell-like spherical 

Velocity gradient present (~ 2 kms _1 pc _1 ) absent 

Velocity dispersion large (~ lkms" 1 ) a along shell-like structure large (~ lkms" 1 )" at central region 

Age of YSO association < 0.5 Myr ~ 1 Myr 



Note. — This table shows a summary of the differences between CWHRs and CWOHRs; spatial distributions, velocity 
structures, and age of YSO association for our target objects. 

a This value is derived from the 2nd moment maps in Figure [71 to II II 



